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ABSTRACT 

We investigate the limits of ground-based astrometry with adaptive optics us- 
ing the core of the Galactic globular cluster M5. Adaptive optics systems provide 
near diffraction-limit imaging with the world's largest telescopes. The substantial 
improvement in both resolution and signal-to-noise ratio enables high-precision 
astrometry from the ground. We describe the dominant systematic errors that 
typically limit ground-based differential astrometry, and enumerate observational 
considerations for mitigating their effects. After implementing these measures, 
we find that the dominant limitation on astrometric performance in this exper- 
iment is caused by tilt anisoplanatism. We then present an optimal estimation 
technique for measuring the position of one star relative to a grid of reference 
stars in the face of this correlated random noise source. Our methodology has 
the advantage of reducing the astrometric errors as ~ l/y/i and faster than 
the square root of the number of reference stars - effectively eliminating noise 
caused by atmospheric tilt to the point that astrometric performance is limited 
by centering accuracy. Using 50 reference stars we demonstrate single-epoch as- 
trometric precision of ~ 1 mas in 1 second, decreasing to <^ 100 /ias in 2 minutes 
of integration time at the Hale 200-inch telescope. We also show that our astrom- 
etry is accurate to ^ 100 /ias for observations separated by 2 months. Finally, 
we discuss the limits and potential of differential astrometry with current and 
next generation large aperture telescopes. At this level of accuracy, numerous as- 
trometric applications become accessible, including planet detection, astrometric 
microlensing signatures, and kinematics of distant Galactic stellar populations. 



1. Introduction 



The benefits of astrometry have long been clear to astronomers. Measurements of par- 
allax and proper motion yield model independent determinations of fundamental quantities 
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like distance and velocity. It is not surprising that astrometry has motivated a wide variety 
of observational programs using many different techniques to answer fundamental questions 
in astrophysics. Potential applications span a wide range of physical scales including: planet 
detection, reconstru ction of the Milky Way's formation, and tests of ACDM cosmology (e.g. 



Unwin et al. 1120071 ). 



The most ubiquitous astr ometric measurement s have been carried out with ground-based 
telescopes in the seeing limit. iMonet et al. I (119921 ) conducted visible light measurements of 
72 stars {V = 15-20) using the 1.55 m US Naval Observatory astrometric reflector. This 
program achieved single epoch measurement precision 4 millia rcseconds (mas), and par - 
allax accuracies ranging from 0.5-3 mas over ~ 5 yr baselines. iPravdo fc Shaklaru (119961 ) 
performed visible light measurements of stars in the cluster NGC 2420 {V = 13-16) and 
achieved single epoch precision of ~ 150 /las in one h our, which motivate d an astrometric 
survey for low-mass companions to nearby stars (e.g. iPravdo et al. II2004I ). More recently, 
200-300 Atas astrometric precision has b een demonstrated with VLT/FORS in the visible 
( iLazorenkd l2006l : iLazorenko et al. 1120071 ). Each of the above programs employed relatively 
narrow-field visible imagers (a few square arcminutes) to perform differential astrometry; 
howev er, the increasing avail ability of wide angle imagers has motivated studies over larger 
fields. [Anderson et al. II2006I performed similar experiments using a 33' x 34' visible camera 
on the ESO 2.2 m telescope, which resulted in 7 mas single-epoch precision. 

Ground-based interferometers provide an alternative method for performing high pre- 
cision astrometry, typically over very narrow fields relative to a single reference star. The 
Palomar Testbed Interferomete r has used phase- referencing to achieve astrometric accura- 



cies «i 100 /ias for a 30^^ binarv (ILane et al. 



2000), and ^ 20 /las over years for binaries with 



separations ^ 1" (IMuterspaugh et al. II2006I ). Due to its 40 cm apertures, this instrument 
is limited to targets with Kg < 6. Large aperture, ground-based interferometers equipped 
with adaptive optics systems, s uch as those at Keck (jColavita &: Wizinowichll2003l ) and the 



[e.g. 



Boden et al. 1120071 ) 



VLT (IGlindemann al. Il2000l ). can perform at similar levels to fainter limiting magnitudes 



Differential astrometric accuracies achieved in both single aperture and interferometric 
ground-based programs are fundamentally limited by atmospheric effects. In the seeing 
limit, single aperture observations suffer from image quality degradation and interferometers 
lose visibility fringe coherence due to atmospheric turbulence. In addition, all ground-based 
programs suffer from systematic effects due to differential atmospheric refraction and optical 
distortions. 



Space-based observatories are one possible method for avoiding the effects of atmo- 
spheric turbulence. Hipparcos was the first space-based mission with astrometric goals. 
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and achieved <; 1 mas astrometry over the mission hfetime on bright targets {V <^ 9 mag; 



Ferryman et al. 1 119971 ) . Currently, the only space-based telescope that can perform high- 
precision astrometry is Hubble. Both the imagers and Fine Guidance Sensor have been 
characterized and well-utilized for astrometry at the <^ 1 mas level (e.g. lAnderson fc King 



2000l . I2OO3I : [Benedict et al. II2003I ). Two complimentary future space missions are aimed 



at achieving levels of astrometric performance 2-3 orders of magnitude below the Hubble 
performance levels. GAIA will catalog roughly one billion stars to ~ 20 mag over the 
entire sky with paralla x accuracies ranging from 10-300 fias depending on the magnitude 
( [Ferryman et al. II2OOII ). The Space Interferometry Mission (SIM) will take a pointed ap- 
proach, and enable microarcsecond (/ias) astrometry on Galactic and extragalactic targets 
dUnwin et allboOTh . 



Ground-based adaptive optics (AO) offer an alternative, more easily accessible, and cost 
effective method for overcoming atmospheric turbulence over small fields (^ arcminute). The 
current generation of astronomical adaptive optics systems provide diffraction limited image 
quality at near-infrared wavelengths. Achieving the telescope's diffraction limit and the 
resulting boost in signal-to-noise ratio prove to be a powerful combination for astrometry. 
These two effects reduce the errors in determining stellar centers, increase the number of 
possible reference stars at small separations, and allow techniques for mitigating systematics 
(e.g. use of narrow-band filters to eliminate chromatic refraction; see 

The marked improvement in wavefront sensor technology and the de velopment of laser 
beacons has rapidly increased the usable sky coverage of these systems (e.g. IWizinowich et al. 
20061 ). The increase in sky coverage, operation in the near- infrared, gain in signal-to-noise ra- 
tio, and the diffraction-limited image quality make astrometry with adaptive optics amenable 
to numerous Galactic applications spanning a wide number of fields: detection of astrometric 
companions, the improved determination of the mass-lurninosit y relation of stars, and the 
formation and evolution of compact objects (lUnwin et al. 1120071 ). 



Here we present an optimal estimation technique appropriate for mitigating the astro- 
metric errors arising in AO observations and demonstrate its potential with multi-epoch 
imaging of the core of the globular cluster M5 using the Hale 200-inch Telescope. We are 
able to achieve <^ 100 //as astrometric precision in 2 minutes, and have maintained this ac- 
curacy over 2 months. In ^we discuss the dominant noise terms that arise in ground-based 
astrometry and the experimental techniques we have adopted to control them. We lay out 
the framework of our reduction model and illustrate its salient properties with a numerical 
simulation in ^ We describe the observations of M5 and the results of applying the optimal 
estimation technique to the data in ^ and ^ This is followed in ^|6] by a discussion of the 
role and potential of adaptive optics in ground-based astrometry with current and future 
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large aperture telescopes. 



2. Astrometric Error Terms in Ground Based Astrometry 

Ground-based optical and infrared imaging observations suffer from a number of errors 
that limit the accuracy and precision of astrometric measurements. Relative to seeing-limited 
observations, the diffraction-limited image quality afforded by adaptive optics modifies the 
relative importance of these error terms. This section describes the four largest effects, and 
indicates observational considerations utilized in this experiment aimed at mitigating them. 



2.1. Differential Tilt Jitter 



With AO, the image motion of the guide star is removed with a flat tip-tilt mirror. 
This stabilizes the image of the guide star with respect to the imager to high accuracy. Any 
residual tip-tilt error is removed in subsequent analysis by calculating only differential offsets 
between the target of astrometry (not necessarily the AO guide star) and the reference stars. 
However, the difference in the tilt component of turbulence along any two lines of sight in 
the field of view causes a correlated, stochastic change in their measured separation, known 
as differential atmospheric tilt jitter. 

More specifically, in propagating through the atmosphere to reach the telescope aper- 
ture, light from the target star and light from a reference star at a finite angular offset traverse 
different columns of atmospheric turbulence that are sheared. Differential atmospheric tilt 
jitter arises from the decorrelation in the tilt component of the wavefront phase aberra- 
tion arising from this shearing effect. This differential tilt leads to a random, achromatic, 
and anisotropic fluctuation in the relative displacement of the two objects. The three term 
approximation to the parallel and perpendicular components of the variance arisi ng from 
differ ential atmospheric tilt jitter, assuming Kolmogorov turbulence, is given by (jSasiela 
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In this equation D is the telescope diameter and Q is the angular separation of the stars. 
The turbulence moments [im are defined as 



sec"^+^e / dhC'^{K)W' 



(2) 
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where h is the altitude, ^ is the zenith angle, and C^{h) is the vertical strength of atmospheric 
turbulence. Typical C^{h) profiles yield cr|| tj !^ 20-30 mas for a 20" binary when observed 
with a 5 m aperture. Note that the variance from differential tilt is a random error, and thus 
is also oc TTj/t, where ttj is the tilt jitter timescale (of order the wind crossing time over 
the aperture; see ^ and t is the integration time. 



2.2. Distortion 

The largest instrumental systematic that limits the accuracy of astrometry in any optical 
system is geometric distortion. These distortions can be stable — resulting from unavoidable 
errors in the shape or placement of optics — or dynamic — resulting from the flexure or 
replacement of optics. 

If geometric distortions are stable, then a number of strategies can be employed to 
mitigate their effect. One method is to mod el the distortion to high accuracy; the most 



notable example is the calibration of HST (e.g. lAnderson fc Kingll2003l ). This is particularly 
important for data sets obtained with multiple instruments or those that use the technique of 
dithering, since knowledge of the distortion is necessary to place stellar positions in a globally 
correct reference frame. Alternatively, one could use a consistent optical prescription from 
epoch-to-epoch by using the same instrument and placing the field at the same location and 
orientation on the detector. Here we use both a distortion solution and a single, consistent 
dither position to achieve accurate astrometry. 

Any changes in the geometric distortion must be tracked through routine, consistent 
calibration. The question of stability is particularly important at the Hale 200-inch, since the 
AO system and the imaging camera (PHARO; see Q are mounted at the Cassegrain focus, 
and PHARO undergo es a few warrn ing /cooling cycles per month (see The PHARO 
distortion solutioij^ by iMetchev ( 2006 ) accounts for changes in the orientation of the telescope 



(which are relatively small for our experimental design), but the overall stability of the system 
is best verified with on-sky data. One of the purposes of the data presented here is to track 
the system stability. We find that the combination of the Hale Telescope, PALAO, and 
PHARO is capable of delivering ^100 /las astrometry. 



see also 



http : / / www . astro. ucla.edu/~ met chev / ao . ht ml 
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2.3. Atmospheric Refraction 



Refraction by the Earth's atmosphere causes an angular deflection of light from a star, 
resulting in an apparent change in its position. The magnitude of this deflection depends 
on the wavelength and the atmospheric column depth encountered by an incoming ray. The 
former effect is chromatic, while the latter is achromatic. The error induced by differential 
chromatic refraction (DCR) has proven to be an im portant, and somet i mes the dominant, 
astrometric limitation in ground-based effor ts (e.g. iMonet et al. Ill992l : iPravdo fc Shakh 
19961 : [Anderson et al. ll2006l : lLazorenkdl2006l ). These studies have shown DCR can contribute 
^ 0.1-1 mas of error depending on the wavelength and strategy of the observations. 

The observations presented here were conducted using a Br-7 fllter at 2.166/im with 
a narrow bandpass of 0.02 /im to suppress differential chromatic refraction. The increased 
signal-to-noise ratio provided by adaptive optics allows sufficient reference stars to be de- 
tected even through such a narrow fllter in a short exposure time. We reach Kg 15 magnitude 
in our 1.4 s exposures through this fllter with the Hale 200-inch (see §3])- In addition, obser- 
vations were acquired over a relatively narrow range of airmass (1.17-1.27) at each epoch to 
minimize the achromatic differential refraction. 

In order to estimate the effect of atmospheric refraction on our data we took the asterism 
in the core of M5 and refracted it to 37 and 32 degrees elevation with the parallactic angles 
appropriate for the obse rvations on 2007 May 2 8 using the slarefro function distributed with 
the STARLINK library (iGubler Sz Tytlerlll998l ). The root-mean-square (RMS) deviation in 
reference star positions between these two zenith angles was <^ 250 yuas and the shift in guide 
star position with respect to the grid (see §1]) was ~ 10 /las. Thus, our consistent zenith 
angle of observations, narrow-band fllter and observations in the near-infrared (where the 
refraction is more benign) make the contribution of this effect negligible for our purposes, 
and we make no effort to correct for it. 



Performing a similar experiment using a the broadband K fllter with a fleld of ~ 5000 K 
reference stars and a ^ 3000 K target would lead to a systematic shift of ~ 100 //as between 
zenith angles separated by 10°, which would be detectable by this experiment. Consequently, 
for observations where broadband fllters are necessary, refraction effects must be considered 
and corrected. 



2.4. Measurement Noise 



In the case of a perfect optical system, a perfect detector and no atmosphere, the astro- 
metric precision is limited one's ability to calculate stellar centers. The centering precision 
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is determined by measurement noise, and we will use the two terms interchangeably. For a 
monopupil telescope the uncertainty is 
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(3) 



(ILindegrenl Il978l ). Adaptive optics allow us to achieve the diffraction limit even in the 
presence of the atmosphere and substantially boosts the SNR over the seeing-limited case 
— thereby decreasing measurement noise and improving astrometric precision. 

In practice, the centering of a given stellar image is limited by spatial and temporal 
variations in the AO point-spread function (PSF). A great deal of time and effort has been 
spent determining the AO PSF and producing software packages to perform PSF fitting (e.g. 



Diolaiti et al. I l2000l : iBrittonI l2006l ). However, any PSF- fitting software package is capable 



of calculating image positions at <^ 0.01 pixel level in a single image. For the observations 
considered here this is ^ 2 mas, a factor of 5 - 10 larger than the measurement noise in 
Equation [3], but it is much smaller than the tilt jitter mentioned in §2.1[ As such, we have 
chosen to use simple and widely available PSF centering software (DAOPHOT; IStetsorull987l : 



see 



3. Observations and Data Reduction 



We observed the the globular cluster M5 on three dates spanning 2 months (see Ta- 
ble [1] for a summary of observations) using t he Hale 200-inch teles cope and the Palomar 
High Angular Resolution Observer fPH ARO: iHavward et al. I l200ll ) assisted by the Palo- 
mar Adaptive Optics System (PALAO; iTroy et al. I l2000l ). The globular cluster M5 was 
chosen for its relatively large distance of ~ 7.5 kpc from the Sun, lo w velocity dispersion of 
K, 5kms~^, and the availability of guide stars near the cluster core (jPryor fc Meylarull993 



Harrislll996l ). This combination of distance and velocity yields an expected cluster dispersion 



Table 1. Observations 
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(mV3) 


(^13/3) 


(„15/3) 


2007-05-28 


05:18:29 - 06:19:29 


890 


1.26 - 1.18 


1.22 


2.34 


l.Ole-5 


3.82e3 


2.77e6 


2007-05-29 


05:58:26 - 06:48:30 


570 


1.19 - 1.17 


1.39 
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2.82e6 


2007-07-22 


03:57:12 - 04:40:36 


630 


1.20 - 1.27 


1.05 


1.66 
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"Calculated at a wavelength of 0.5 ^m. These quantities scale as A^/^ and A^/^, respectively. 
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Fig. 1. — Top left: Image of the core of the globular cluster M5 in 1.4 seconds through the 
narrow-band Br-7 filter. The AO guide star is labeled as star 'A', and is one of 82 detected 
stars in the image. The additional plots show the measured x-y angular separation of each 
pair of stars denoted by the arrows in 600 images taken on 28 May 2007. These plots show the 
clear signature of anisotropic differential atmospheric tilt jitter as predicted from Equation [1] 
The measured (red) and predicted (green) 1 a error ellipses are over-plotted. We see that 
temporal averaging over the 1.4 second exposure time has reduced the measured variance 
with respect to that predicted from the DIMM/MASS measurements and Equation [H 
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of only 140 yuas yr~^, or 20 yuas over our 2 month observing span. We acquired 400-600 images 
per night. A typical image can be found in Figure [TJ The guide star is a red giant branch 
member of the globular cluster with V ~ 12.6 magnitude fISandquist Sz Boltdl2004l ). The 
cluster was imaged through the narrow-band Br-7 filter (central wavelength is 2.166 fim and 
bandpass is 0.02 /zm) using the 25" x 25" narrow-field channel (0.025" pixel"^), which over 
samples the 87 mas diffraction-limited PSF. The brightest star filled the detector to 10% of 
the maximum well-depth in the 1.4 sec exposure time — well within the linear regime of the 
detector. 

Contemporaneous measurements of the atmospheric turbulence profile were acquired 
with a differential image motion monitor (DIMM) and multi-aperture scintillation sensor 
(MASS), which have been deployed as a single unit in a dome at Palomar Observatory 



(IThomsen. Britton fc Pickled 120071 : iKornilov et al. 1 120071 ) . These turbulence profile mea- 



surements permitted an independent estimate of the magnitude of differential tilt jitter 
(computed using Equation [1]) . 

We processed the raw images by subtracting dark frames and removing bad pixels from 
the analysis. Flat-field calibration was performed using twilight sky flats. Sky subtraction 
was accomplished by forming the median of the dithered frames taken outside of the cluster 
and subtracting this median from each exposure. The photometry and astrometry of each 
star was extracted using PSF-fitting as implemented by the DAOPHOT package in PyRAF@. 
DAOPHOT is not optimized for astrometry (see e.g. Anderson fc King 200ol ). and since 
our measurement model reduces the noise due to atmospheric turbulence, our single epoch 
precision could be improved with a more careful centering technique (see §4.41 and §5.2p . 
However, our astrometric accuracy over 2 months is not limited by this choice (see §5.31) . We 
used the 4 brightest stars in the field to derive a model PSF that is assumed to be constant 
over the field, and calibrated the image zeropoints using 2MASS and find that the guide star 
has Kg K, 9.1 magnitude. 



4. Grid Astrometry for Ground-based Adaptive Optics Observations 

After controlling for distortion and atmospheric refraction, the dominant astrometric 
errors are caused by differential atmospheric tilt jitter and measurement noise. In this sec- 
tion we present a general framework for measuring the position of a star relative to a grid 
of reference stars in the face of these noise sources. This framework has two key ingredients. 
The first is the covariance matrix (Sd), which encapsulates the relevant statistical uncer- 



^PyRAF is a product of Space Telescope Science Institute, which is operated by AURA for NASA. 
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tainties for astrometry with adaptive optics. The second is the weight matrix (W), which 
determines how the differential measurements between the target star and the reference stars 
are combined to calculate the target's position relative to the grid. 



4.1. Measurement Model 



The fundamental quantity in differential astrometry is the measured angular offset be- 
tween a pair of stars. We will denote the angular distance between two stars, i and j, as 
dij. Since dij is measured from an image, we will denote its components in the Cartesian 
coordinate system of the detector, simply 



di 



where we have introduced the notation Xij 



Xj 
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(4) 



Xi and likewise for y. The variance in the 



the angular separation between two stars is given by 
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(5) 



where cxy is the sum of the squares of the centering errors of each star parallel to the axis 
connecting the pair (and similarly for the perpendicular variance), and the remaining terms 
are as defined in §2.11 

Measurement of the offset between the target star (which we will denote with a subscript 
i = 0) and each of the N reference stars results in a set of vectors, doi. For simplicity, we 
will write these measured offsets as a single column vector. 



[Xoi, ■ ■ ■ ,XoN,yoi, ■ ■ 



(6) 



The goal of differential astrometry is to use d to determine the position of the target star 
with respect to the reference grid of stars at each epoch. 

There are many possible ways to construct the position of the astrometric target from 
a given d. Here we use the most general linear combination of the angular offsets, namely 



P 



Wd, 



(7) 



where W is the 2 x2A^ weight matrix, given by 



W 
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We have used the notation Wxy^i to denote the weighting of the offset from the target star 
to star i in the y direction used to determine the x component of the target's position, p. 
For example, for a standard average of the x and y measurements to calculate p, we would 
assign all the w^^fii = Wyy^i = 1/N and w^yfii = Wy^fii = 0. 

In principle, we are free to assign weights in any manner we please. However, we find it 
convenient to choose the weights such that they satisfy 

Wxxfii = 1 ' X] ^yyfii = 1 ' X] ^^yfii = ' X] ^v^fii = (9) 

i i i i 

These constraints ensure that the components of p have physical units (e.g. pixels or arcsec- 
onds) and that its components are measured in the same coordinate system as d (presumably 
the detector coordinates). As a consequence, p represents the position of the target star in 
the sense that a proper motion of the target, e, with respect to the fixed grid between two 
epochs will cause a change, p ^ p + e. 

In order to determine if any change in p over time is meaningful we must understand 
its statistical properties. Both differential tilt jitter and measurement errors are assumed to 
follow Gaussian statistics, so that each instance of target-reference grid offset measurements, 
d, is drawn from a multivariate normal probability distribution: 

where is the covariance matrix, and the bars above symbols denote using the average 
value of each matrix entry. 

The statistics of p follow in a straightforward manner from Equation [TO] given our choice 
in Equation [71 Since p is a linear function of d, each p is also drawn from a multivariate 
normal probability distribution with covariance matrix 

Sp = W^SdW, (11) 

and the uncertainties of p are described by the eigenvectors and eigenvalues of Sp. Thus, 
our goal of optimally determining the target's position requires calculating the covariance 
matrix, Sd, from data or theory, and choosing W to minimize the eigenvalues of Sp. 
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4.2. The Covariance Matrix 



We have chosen to measure positions and offsets in the Cartesian coordinates of the 
detector, so the form of the covariance matrix, given our definitions above, is 

/ ((Axoi)2) ••• ((Axoi)(Axo7v)) ((Axoi)(Ayoi)) ••• ((A2;oi)(Ayo7v)) \ 



((Axo7v)(Ayoi)) 



y symmetric 
where we have written Ax,- 



X 



((AxoAf)(AyoAr)) 
((Ayoi)(Ayoiv)) 

{{AyoN?) J 
to simplify the notation (likewise for y). 



(12) 



The total covariance matrix has contributions from centering errors and differential 
atmospheric tilt jitter. Since these contributions are independent, the total covariance matrix 



can be written 



+ Stj, and each term can be derived separately. 



4-2.1. The Covariance Matrix for Measurement Noise 

In the absence of differential tilt jitter it is straightforward to construct the covariance 
matrix for measurement noise alone, Smeas- The diagonal terms can be written 

= <o^ = <o + (13) 

where ax^i and ax,o are the the uncertainties in determining the x-position of star i and the 
target star, respectively. For the off-diagonal terms (AxoiAxoj) we can use the fact that 

(AxcAxo.) = l({A4 + Ax^,-[Axo.-Axo,f}) 
= ^{(Ax^.) + (A4)-(Ax^,)} 

\( 2 _|_2 _2n 

= <o- (14) 

where we have used only algebra and the definitions above. Equation [H] is the obvious 
result of the fact that the measurements of the target star's coordinates are common to all 
differential measurements, and so its uncertainty appears in all the off-diagonal covariance 
terms, (AxojAxoj) and (AyoiAyoj)- However, the cross-terms involving both x and y (e.g. 
(AxoiAyQj)) vanish because a^^ and a^^o uncorrelated for measurement noise alone. 
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4-2.2. The Covariance Matrix for Differential Tilt Jitter 

The covariance matrix for differential atmospheric tilt jitter between a pair of stars is 
diagonal when written in an orthogonal coordinate system with one axis lying along the 
separation axis of the binary. We see from Equation [1] that it can be written as 

where d\\ and d±_ are the angular offsets parallel and perpendicular to the axis connecting 
the pair of stars, respectively. 

For a general field of stars, no coordinate system exists that diagonalizes the full tilt 
jitter covariance matrix, Stj. But, we can begin computing the entries by rotating Spair 
into our x-y coordinates via R^SpairR, where 

R=f """"t '"tV (16) 
Y — sm (j) cos (j) J 

The result is 

rT^ ^ _ ( ^|oi cos^ 0oi + ^i,oi sin^ 0oi (^|oi ' ^i,oi) cos 0oi sin 0oA . . 
\ (^i|,oi ~ ^l,oJ COS 0oi sm 0oj o-pi 0oi + oi cos^ 0oi y 

where (/)oj is the angle between cioi and our arbitrary Cartesian system measured counter- 
clockwise from the a;-axis, and we have introduced the notation that the uncertainty parallel 
to dij is (T||^jj and uncertainty orthogonal to di^ is o^^i^ as calculated from Equation [TJ Thus, 
we can identify the diagonal terms 

^^^lil = ^ifli cos^ 0oi + (^l,oi sin^ <Poi, 

i^Voi) = (^loi sin^ 0oi + orlfli cos^ 0oi, (18) 

and 

{AxoiAyoi) = (crj - cr^ cos 0oi sin 0oi- (19) 
For the off-diagonal terms {AxoiAxoj) we notice that (as used in Equation [1] 

(A.,.A.,,) = 1({A4 + A4-|A.„.-A.„,in> 

1 

2 
1 



{(A4) + - ([Axo. - Axo.f )} 



2 l^y,oi cos" (poi + al oi sin" (poi + a^ Qj cos^ + al ^j sin^ 
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where in the last step we have used the fact that Xij = xoi — xqj and the relations in 
Equations [T8] and [T9l The quantities (A|/oiAyoj) can be obtained by interchanging sine and 
cosine in Equation [2U1 

For the remaining off-diagonal terms (AxoiAyoj) we can use the fact that 

(AxoiAyoj) = {Axoi[Ayoi + Ayij]) 

= (AxoiAyoi) + {AxoiAyij]) 

= {AxoiAyoi) + {[Axoj - Axij]Ayij) 

= (AxoiAyoi) + {AxojAyij) - {AxijAyij) 

= (AxoiAyoi) + {Axoj[Ayoj - Ayoi]) - {AxijAyij) 

= {AxoiAyoi) + {AxojAyoj) - (AxijAyij) - (AxojAyoi). (21) 

Rearranging gives 

{AxoiAyoj) + (AxojAyoi) = (AxoiAyoi) + (AxojAyoj) - (AxijAyij). (22) 

All the terms in the right-hand side are known from Equation [T9|, and further investigation 
shows that the two terms on the left-hand side are equal. So, Equations [131 [HI [HI [IHl [201 
and [221 contain all the information required to construct the full covariance matrix, Sj. 



4.3. The Optimal Weight Matrix 

The optimal choice of weights in Equation [3] are those which minimize the eigenvalues 
in Equation [Til For a 2 x 2 symmetric matrix the sum of the eigenvalues is the trace of the 
matrix, so our problem becomes one of minimizing the trace of Sp subject to the constra ints 
in Equation [HI Specifically, we will use the method of Lagrange multipliers (lBettslll980l ) to 
find the optimal weights, W, that minimize the quadratic equation 



Tr(Sp) = iw^SW, 



(23) 



where 



is a AN X AN matrix, and 



Sd 
Ed 



(24) 



W — [Wxx,Ol! ■ ■ ■ ! Ul^xfiN, Ulxyfil, ■ ■ ■ , WxyfiN, Wy^fil, " " " , Wy^fiN, Ulyyfily ' ' ' ' '^vyfiM (25) 
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is a vector of length 4N. Note that W has identical entries as W in Equation [8|, it is just 
written as a single vector to cast the minimization problem into a single equation fl23|) . We 
want to find the extrema of Equation [22] subject to the linear constraints in Equation [HI 
which can be written 

CW' = V, 



(26) 



where we define the 4 x 4N symmetric matrix 

1 ■ 



C 



1 

1 



sym 



and 



V 





1 

1 



1 




1 







1 

1 







1 



(27) 



(2^ 



In this framework the optimal weights are those that solve the system of linear equations 

(29) 
can be 



" s 






" W " 




" " 


c 







A 




V 



Here A are the Lagrange multipliers, which will not be used further. Equation 
solved via a matrix inversion. 

Note that the general constraints on the weights we have written in Equation M and 
have two somewhat unintuitive features. The first is that the y measurements are sometimes 
used to compute the x position and vice versa. The other property is that they allow for 
negative weights, meaning that in some cases certain measurements will be subtracted in 
calculating the position of the astrometric target, p. These two facts conspire to exploit 
the natural correlations inherent in the data. The fiexible and possibly negative weights 
essentially allow the reference grid to be symmetrized, thereby using the known correlations 
to cancel noise so as to minimize the variance in p. 



4.4. Numerical Simulations 



As indicated in the above analysis, the single epoch uncertainty in the location, p, of the 
target relative to the grid of reference stars is represented by the eigenvalues and eigenvectors 
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of the 2x2 matrix Sp (Equation [TT]) . This matrix itself depends on the distribution of 
reference stars, the precision of centering measurements, and the degree of noise correlation 
due to differential tilt through the matrix Sj- In this way, the intrinsic precision of the 
measured value of p depends on these three factors. 

To ascertain the behavior of Sp with the density of available reference stars, we per- 
formed a series of numerical simulations. In each simulation, N {2 < N < 100) stars were 
randomly distributed throughout a 25" x 25" field of view. We assumed the target was a 
bright star in the middle of the field with centering error of 0.5 mas and the reference stars 
were fainter, drawn from a Gaussian distribution with mean centering error of 2 mas and 
a standard deviation of 1 mas (somewhat analogous to the situation for the guide star in 
M5; see ^S])- The full covariance matrix, Sj? was computed for each stellar configuration 
assuming these centering errors, the typical turbulence profile above Palomar Observatory, 
and a 1.4 sec exposure time. 

In the first simulation, Sd was contracted as in Equation [11] using standard averaging 
for W {wxxfli = '^yyfli = V^i ui^yfii = Wyxfii = 0). For the second simulation, Sd was 
contracted using the optimal W as calculated using the prescription in §4.3[ In each case, 
the geometric mean of the two eigenvalues of the resulting matrix, Sp, were computed to form 
an estimate of the single epoch measurement precision of p. To average away random effects 
arising from the particular geometry of the random distribution of stars, each numerical 
simulation was repeated for 100 random distributions of stars for each value of N, and these 
were averaged to generate a mean value for the single epoch measurement precision. 

The resulting values for the single epoch measurement precision of p are shown in 
Figure [2] as a function of the number of reference stars, along with the contributions of mea- 
surement noise and differential tilt jitter. In both simulations the error due to measurement 
noise decreases as A^"'''^. However, in the limit of an infinite number of reference stars, this 
error asymptotes to the target star's measurement error. The rate at which the measure- 
ment noise decreases to this value depends on the distribution of reference star measurement 
errors. 

The important distinction between the two simulations is the contribution of tilt jitter 
to astrometric performance. In the simulation utilizing standard averaging, there is very 
little gain with increased stellar density {N~^'^^), and tilt jitter dominates the error budget. 
However, the optimal estimation algorithm rapidly (A^~° '') eliminates the contribution of 
differential tilt by taking advantage of the correlations inherent in Sd and the flexibility to 
symmetrize the reference field through the choice of weights. 
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Fig. 2. — Top: Simulated astrometric precision as a function of the number of reference stars 
using standard averaging (solid line). The total astrometric precision has contributions from 
the measurement noise (dash-dotted) and the differential atmospheric tilt jitter (dashed line). 
Here the measurement noise was taken to be 0.5 mas for the AO guide star, and the reference 
stars were drawn from a normal distribution with a mean of 2 mas and standard deviation of 
Imas. The tilt jitter is that expected in a 1.4 sec exposure assuming the turbulence profile 
measured on the night of 2007 May 28 at Palomar Observatory (see Table [1]). Bottom: 
Simulated astrometric precision as a function of the number of reference stars using optimal 
weighting (lines as above). By using optimal weighting based on the covariance matrix, the 
effect of atmospheric noise is reduced to values less than measurement noise. 
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5. Analysis and Results 

In the analysis that follows we will use the measurement model described in §11 For a 
given target star, we will calculate the differential offsets with respect to the grid stars to 
generate a value of d for each image at each epoch (Equation [6]). We then use either these 
data or the theory in §4.21 to generate the full covariance matrix, S^. From S^, we use 
the prescription in §4.31 to calculate the the optimal weights, W. These weights are used 
to combine the differential offsets to generate the target star's position, in each image 
via Equation [71 The statistics of the positional measurements are then described by the 
covariance matrix, Sp, from Equation [TTl 



5.1. Differential Tilt Jitter 

In order to test our expectation that tilt jitter dominates the astrometric error, we 
calculate the RMS of the angular offsets for pairs of stars in the field (Figures [H and [3]). 
These results clearly show the characteristic signature of differential tilt. Namely, the RMS 
separation along the axis connecting the two stars is larger than that of the perpendicular 
axis by a factor of ^ ^/3. However, the magnitude of the tilt jitter is smaller than the 
theoretical expectations, which suggests that some of the tilt jitter has been averaged away 
in the 1.4 sec exposure time. 

We have no direct measurement of the wind speed profile over the telescope to calculate 
the expected tilt jitter timescale. Instead, we fit the observed afj and angular offsets using 
the model in Equation [5l with t = 1.4 seconds. The best fit values are CTmeas.ij ~ 2 mas and 
t/rxj ~ 7. This implies that the characteristic timescale for tilt jitter is ~ 0.2 seconds, 
resulting in a wind crossing time of 25msec~^. Turbulence at higher altitudes contributes 
most to the differential atmospheric tilt j itter, and this velocity is typical of wind speeds in 



the upper atmosphere (iGreenwoodl 119771 ) . It is also clear from the figure that a number of 
stars have measurement noise that is much less than 2 mas, thus this number should only be 
taken as characteristic of the faint stars. 



5.2. Astrometric Precision 

The astrometric precision achieved in a single epoch is an important diagnostic of the 
measurement model. On a given night for a given star, we investigate the use of both the 
^ 500 images and the theory in §4.21 to calculate Sd, leading to the optimal weights. We 
then apply this weight matrix to the measured offsets to compute the target's position in 
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Fig. 3. — Top: RMS deviation in the distance between pairs of stars in the direction parallel 
to their separation axis on 2007 May 28 (filled circles). The jitter predicted (assuming no 
measurement noise) from the measured turbulence profiles and Equation [1] (dashed-line) is 
far larger than the measured jitter, indicating that some tilt jitter has been averaged away 
in 1.4 sec. The best fit model (Equation [5]) including averaged tilt jitter and measurement 
noise indicates that the tilt timescale is ~ 0.2 sec (solid line). Bottom: As above, but in this 
case the separations and predictions are for the direction perpendicular to the separation 
axis. The expected RMS for the perpendicular direction is lower by the expected factor as 
seen in Equation [H Note that not all pairs include the AO guide star. 
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Fig. 4. — Allan deviation in the guide star position as a function of integration time for 2007 
May 28 (solid curve), May 29 (dash-dotted) and July 22 (dashed). The astrometric precision 
scales as t-o-5i±o.o8^ ^^^^ ^^le covariance matrix and optimal weights were derived from data. 
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each image — resulting in a timeseries in each component of p for each epoch. The properties 
of each timeseries are best explored by computing its Allan deviation (also known as the 
square root of the two-sample variance). The Allan deviation is calculated by dividing a 
timeseries into chunks, averaging each segment, and computing the RMS of the resulting, 
shorter timeseries. If the timeseries is dominated by random errors, its Allan deviation will 
decrease as 1/ y/tavg, where tavg is the length of each chunk. It is also necessary to have 
sufficiently many segments so that an RMS calculation is meaningful. Here the longest 
timescale probed is ~ 2-3 minutes for each 10-15 minute timeseries. 

We compute the geometric mean of the Allan deviation in each dimension as a function 
of the averaging time for the AO guide star in Figure H] after computing the covariance 
matrix from data. After 1.4 seconds the guide star's positional precision is ~ 600 /las. The 
precision subsequently improves as /ias after 2 minutes, and has yet to hit 

a systematic floor. This suggests a precision of ~ 30 /las for the full 10-15 minutes data set, 
assuming that no systematic limit is reached in the interim. 

This level of precision is not limited to the AO guide star; similar performance is obtained 
on other stars in the core of M5. In Figure we show the astrometric precision obtained 
on 2007 May 29 after 2 minutes for all detected stars as a function of their Kg magnitude. 
Precision below 100 /ias is achieved on targets as faint as Kg ~ 13 magnitude using a narrow- 
band filter and 1.4 sec individual exposures. This demonstrates the substantial signal-to-noise 
ratio benefit afforded by adaptive optics. 

The astrometric precision shown in Figure [5] resulting from the theoretically determined 
covariance matrix and optimal weights is ~ 300 /ias after 2 minutes for stars with Kg <s 13 
magnitude. This level of precision is substantially better than the performance of simpler 
weighting schemes, but it is a factor of 2-4 worse than using the data to calculate the 
covariance matrix and weighting. There are several possible reasons for this reduction in 
precision. The first is that we have only used estimates of the measurement noise for each star 
used to calculate Smeas- Secondly, the turbulence profile used to construct Stj is estimated 
from the average C^{h) seen by the DIMM/MASS. This unit is located 300 m from the Hale 
telescope and uses Polaris to estimate the turbulence profile. As a consequence, there could 
be important differences between the measured atmospheric turbulence and that encountered 
by the light from M5. Finally, we have not attempted to capture the time variability of the 
turbulence, having used only the average values. 

In Figure we investigate the improvement of the AO guide star astrometry with the 
number of reference stars. We drew random subsets of the available grid stars, computed 
Sd from the data, calculated the optimal weights, and show the geometric mean of the 
eigenvalues of Sp. To average over the geometry of a particular draw, we repeated this 
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process 10 times for each value of and averaged the results. We see that the precision 
rapidly decreases as ]\f-o-^o^o-03 _ xhis is slightly faster than our simulations predict for 1.4 sec 
of integration time. However, as noted above, our simulations are meant to approximate M5, 
but do not capture the true distribution of stellar measurement errors (which are difficult to 
decouple from tilt jitter) or any evolution in atmospheric turbulence during the observation. 



5.3. Astrometric Accuracy 

The goal of astrometry is to measure the position of the target star over many epochs. 
Astrometrically interesting timescales range from hours to years. Clearly, the optical systems 
must be stable over these spans for astrometry with AO to be viable. There are several obsta- 
cles that could render the single-epoch precision obtained in §5.21 meaningless. For example, 
PHARO is mounted at the Cassegrain focus which results in flexure of the instrument as 
the telescope tracks, and undergoes warming and cooling cycles between observing periods 
(typically twice per month) that could cause small changes in the powered optics. Either of 
these facts could alter the geometric distortion, and make astrometric measurements unre- 
peatable. In order to probe the system stability, we have designed our experiment to be as 
consistent as possible, and it has spanned many removal and reinstallations of PHARO over 
2 months. 

In order to investigate the accuracy of the M5 measurements we flrst measured and 
corrected the small rotational (^ 0.04°) and plate scale (^ 10"^) changes between the the 
May 29 and July 22 data and the May 28 images. We also calculate the optimal weights for 
a given star on all three nights, and average them to create one weighting matrix to use for 
each epoch. This is not strictly optimal, since each night has different turbulence conditions 
for example, but it ensures that the scenario that p ^ p + e. We see in Figure [7] that the 
measured position of the AO guide star is accurate from epoch-to-epoch at the ~ 100 //as. 
The error ellipses are those estimated by continuing to extrapolate the precision found in 
Figure m by l/-\/t to the full 10-15 minute timeseries. This is an impressive level of accuracy, 
but unfortunately is a factor of 3 worse than our expectation. It suggests that there is some 
instability, likely in the distortion, over the two months that limit the astrometric accuracy. 

The other stars in M5 show a similar level of astrometric accuracy (Figure [H]) up to 
Ks ~ 13 mag. This limit can certainly be pushed considerably fainter with in increased 
integration time or a larger aperture. The achievement of such high levels of astrometric 
performance on faint targets, given the modest time investment, short integration time and 
narrow-band fllters, illustrates the substantial signal-to-noise ratio gain and potential for 
astrometry enabled by AO. 
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Fig. 5. — The astrometric precision (Allan deviation after 2 minutes) using the theoretical 
covariance matrix (open diamonds) and the covariance matrix from data (filled circles) as 
function of Kg magnitude for all 82 detected stars on 2007 May 29. The precision in both 
cases is essentially constant for Ks ^ 13 mag. However, the astrometric precision for the 
theoretical Xld is a factor of 2-4 times larger than when calculated from data. 
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Fig. 6. — Allan deviation in 1.4 seconds for the AO guide star's position as a function of the 
number of reference stars on 2007 May 28 (solid curve), May 29 (dash-dotted), and July 22 
(dashed). The astrometric precision scales as jv~'^'^°^'^'°^. 
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Fig. 7. — The position of the AO guide star in an arbitrary coordinate system on three dates: 
2007 May 28 (sohd), 2007 May 29 (dashed), and 2007 July 22 (dash-dotted; see §53]). The 
error circles are inferred by averaging the covariance matrix measured from the data and 
extrapolating to the total 10-15 minute integration time as 1/v^ (e.g. see Figure Hj). The 
positions agree at the ^ 100 yuas level — a factor of 2-3 larger than the expected dispersion. 
This discrepancy indicates that some systematic errors have occurred between epochs, most 
likely optical distortion. 
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Fig. 8. — Astrometric accuracy (geometric mean of the RMS in each coordinate over the 
three epochs) versus Kg magnitude. The level of accuracy is <^ 100 yuas, and is essentially 
unchanged for Kg ^ 13 mag. However, this is a factor of 2-3 above our expectations from 
the achieved astrometric precision, suggesting a systematic limitation between epochs. 
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6. Discussion and Conclusions 

Here we have presented a technique for performing high-precision grid astrometry us- 
ing ground-based telescopes equipped with adaptive optics systems. With this technique, 
the effects of distortion and atmospheric dispersion that give rise to systematic errors are 
mitigated through the design of the experiment. Random errors arising from differential 
tilt jitter and measurement noise are minimized through the use of an optimal estimation 
scheme that accounts for the correlated noise statistics through the covariance matrix Sd. 
The experimental results obtained on the Hale 200-inch Telescope have demonstrated single 
epoch astrometric precision of <^ 100 //as in 2 minutes and multi-epoch astrometric accuracy 
at the same level. This level of precision is comparable to that afforded by ground-based 
interferometry, and is better than the precision obtained in seeing-limited programs on single 
apertures. 

The simulation of astrometric precision afforded by the optimal weighting scheme shown 
in Figure [2] illustrates that measurement noise is the dominant residual astrometric error on 
a 5 meter telescope for stellar fields that contain more than a few reference stars. The scaling 
laws for differential tilt jitter and measurement noise (-D~^) indicate that on larger 

aperture telescopes measurement noise will represent a smaller fraction of this residual error. 
This effect is illustrated in Figure [HI which shows the RMS error between pairs of stars for 
a range of telescope apertures and angular separations. 

The values in Figure [9] assume that tilt jitter is resolved by sufficiently short exposures. 
Longer exposure times will certainly reduce the differential tilt jitter by 1/Vi, but the mea- 
surement noise will also be decreased by this factor (for a given stellar brightness). The 
implication being that if tilt jitter dominates for short exposure times, it will continue to 
dominate longer exposures. 

In situations where fewer images are available, either due to time constraints or longer 
exposure times per frame, it is difficult or impossible to effectively calculate the covariance 
matrix directly from the data. Our results show that independent measurements of the 
turbulence profile, for example from a DIMM/MASS unit, are sufficient to calculate Sd, 
and result in astrometric precision within a factor of 2-4 of the levels achieved using the 
data itself. Thus, the astrometric applications of turbulence sensors are two-fold; t hey can 



be us ed to independently assess astrometric data quality, and predict the AO PSF (IBritton 



20061 ). 



The scaling laws presented throughout this paper indicate a substantially improved 
astrometric performance on large aperture telescopes equipped with adaptive optics. We 
have used the measured performance on M5 with the Hale Telescope combined with these 
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Fig. 9. — The RMS separation between a pair of stars versus angular offset and aperture 
diameter. We have assumed the turbulence profile on 2007 May 28, a measurement error of 
1 / a/2 mas for each star for a 5 m telescope, and included the geometric mean of each com- 
ponent of Equation [H Relative to Palomar (solid curve), there are substantial astrometric 
gains to be made for larger 10 m (dash dotted) and 30 m (dashed) telescopes due to the 
reduction of both measurement noise (the y-intercept; oc D~'^) and tilt jitter (oc D""^/^). 
Because measurement noise falls off more quickly with D, tilt jitter becomes the dominant 
source of astrometric error for large aperture telescopes. 
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Fig. 10. — Astrometric precision as a function of aperture size and stellar density. We have 
used Equation [301 with the assumptions of the Palomar turbulence profile, a 25" x 25" field 
of view, the M5 brightness distribution, and photon noise limit as described in ^ The 
astrometric precision demonstrates a very favorable scaling law with aperture diameter, and 
suggests orders of magnitude improvement in precision may be available using large aperture, 
AO equipped telescopes. In practice, the level of astrometric accuracy will depend on the 
extent to which current and future facilities can characterize and control systematic errors. 
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scaling laws to predict the astrometric performance of a single conjugate AO system as 
a function of aperture diameter and number of reference stars. The relationship can be 
summarized using the results of simulation and data analysis as 



J2 _ 2 I 2 
'^tot ~ '^mcas'T'^TJ 



1.4 sec' 

V 



2 mas 



\ 0.3 /-. \ 2 

1 \ /5m 



N \ D 



2 



2 mas 



o \ 0.7 /p. N 7/6' 

2 \ / 5 m \ 



N J \ D J 



2 



(30) 

This equation assumes that measurement error is dominated by photon noise (oc D~^) and 
the other dependencies (field of view, stellar brightness distribution, turbulence profile) are 
identical to those for the M5 experiment. 

Figure 10 shows the resulting estimates for astrometric precision as a function of aper- 
ture diameter and number of reference stars for a 2 minute exposure. These predictions 
demonstrate that limits to astrometric precision arising from random errors (dominated by 
tilt jitter) lie below 10 yuas for 30 m telescopes. However, very careful characterization and 
control of systematic errors will be be required to achieve this level of precision in an actual 
experiment. The extent to which systematic errors can be eliminated will distinguish the 
scientific goals that can be accomplished with ground-based facilities from those that require 
a space-based solution. 

Facilities: Hale (PALAO/PHARO) 
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